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DISCRETE APPROXIMATIONS TO REFLECTED 
BROWNIAN MOTION^ 

By Krzysztof Burdzy and Zhen-Qing Chen 

University of Washington 

In this paper we investigate three discrete or semi-discrete ap- 
proximation schemes for reflected Brownian motion on bounded Eu- 
clidean domains. For a class of bounded domains D in R" that in- 
cludes all bounded Lipschitz domains and the von Koch snowflake 
domain, we show that the laws of both discrete and continuous time 
simple random walks on Dn2^'°Z" moving at the rate 2^^*° with sta- 
tionary initial distribution converge weakly in the space D([0, 
equipped with the Skorokhod topology, to the law of the stationary 
reflected Brownian motion on D. We further show that the following 
"myopic conditioning" algorithm generates, in the limit, a reflected 
Brownian motion on any bounded domain D. For every integer fc > 1, 
let {X^2-fc li = 0, 1, 2, . . .} be a discrete time Markov chain with one- 
step transition probabilities being the same as those for the Brownian 
motion in D conditioned not to exit D before time 2~*. We prove 
that the laws of X* converge to that of the reflected Brownian mo- 
tion on D. These approximation schemes give not only new ways 
of constructing reflected Brownian motion but also implementable 
algorithms to simulate reflected Brownian motion. 

1. Introduction. Let n >1 and D C M" be a domain (connected open 
set) with compact closure. Consider a reflected Brownian motion (RBM in 
abbreviation) Y in D. Heuristically, RBM in D is a continuous Markov 
process Y taking values in D that behaves like a Brownian motion in M" 
when Yt £ D and is instantaneously pushed back along the inward normal 
direction when Yj € dD. RBM on smooth domains can be constructed in 
various ways, including the deterministic Skorokhod problem method, mar- 
tingale problem method, or as a solution to a stochastic differential equation 
with reflecting boundary conditions (see the Introduction of [5]). When D is 
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nonsmooth, all the methods mentioned above cease to work. On nonsmooth 
domains, the most effective way to construct RBM is to use the Dirichlet 
form method. The RBM constructed using a Dirichlet form coincides with 
RBM constructed using any other standard method in every smooth domain. 

It is natural to try to construct RBM in a nonsmooth domain using a 
sequence of approximations that can be easily constructed themselves. One 
such approximating scheme was studied in [2] and [5], where the processes 
approximating RBM in a nonsmooth domain D were RBMs in smooth do- 
mains increasing to D. In this paper we will consider processes approximat- 
ing RBM in D that are defined on the same state space or a discrete 
subspace of D. We consider the subject interesting and important in itself, 
but we also have a more concrete and direct motivation — it comes from 
two recent papers on multi-particle systems [3, 4]. In each of these papers a 
large population of particles is trapped in a domain. All particles perform 
independent random walks or Brownian motions and an appropriate mech- 
anism keeps all particles inside the domain, which is assumed to be regular 
in a certain sense (see [3, 4] for more details). The intuitive interpretation 
of those results and their extension to less regular domains would require 
an invariance principle for the reflected random walk and similar processes 
in domains with nonsmooth boundaries. As far as we can tell, such results 
are not available in literature and they do not follow easily from published 
theorems. 

In this paper we investigate three discrete or semi-discrete approximation 
schemes for reflected Brownian motion. The first two approximations involve 
random walks and we prove that they converge to reflected Brownian motion 
in a class of bounded nonsmooth domains in M" that includes all bounded 
Lipschitz domains and the von Koch snowflake domain. The third scheme is 
based on "myopic" conditioning and it converges to the reflected Brownian 
motion in all bounded domains. We will now describe these schemes in more 
detail. 

Let D be a bounded domain in whose boundary dD has zero Lebesgue 
measure. Without loss of generality, we may assume that G -D. Let be 
the connected component of H 2~'^Z"' that contains with edge struc- 
ture inherited from I^^TU^ (see the next section for a precise deflnition). 
We will use Vkix) to denote the degree of a vertex x in L'fc. Let and 
be the discrete and continuous time simple random walks on mov- 
ing at the rate with stationary initial distribution m^, respectively, 
where mk{x) = ^^2'^^"^ . We show that the laws of both {X^,k > 1} and 
{Y^, A; > 1} are tight in the Skorokhod space D([0, oo), M") of right continu- 
ous functions having left limits. We show (see Theorems 2.4 and 3.3 below) 
that if D satisfles an additional condition (1.1) below, which is satisfled by 
all bounded Lipschitz domains and all bounded uniform domains (see below 
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for the definition), then both {X'', A; > 1} and {Y'', k> 1} converge weakly 
to the stationary reflected Brownian motion on D in the Skorokhod space 
D([0,1],M"). 

The last of our main theorems is concerned with "myopic conditioning." 
We say that a Markov process is conditioned in a myopic way if it is condi- 
tioned not to hit the boundary for a very short period of time, say, units 
of time, where k is large. At the end of this period of time, we restart the 
process at its current position and condition it to avoid the boundary for an- 
other period of 2"^^ units of time. We repeat the conditioning step over and 
over again. Intuition suggests that, when is very small and the process is 
far from the boundary, the effect of conditioning is negligible. On the other 
hand, one expects that when the process is very close to the boundary, the 
effect of conditioning is a strong repulsion from the boundary. These two 
heuristic remarks suggest that, for small 2^^, the effect of myopic condi- 
tioning is similar to that of reflection. A more precise description of myopic 
conditioning of Brownian motion is the following. For every integer A; > 1, 
let {Zj^^f,, j = 0,1,2, .. .} be a discrete time Markov chain with one-step 
transition probabilities being the same as those for the Brownian motion 
in D conditioned not to exit D before time 2~^. The process can be 
defined for t G [(j — 1)2"'^, j2~'^] either as the conditional Brownian motion 
going from to without leaving the domain D or as a linear 

interpolation between Z^^^-^s^^^^ and Z^^_k- We prove in Theorems 5.1 and 

5.6 below that, for any bounded domain D, the laws of Z^ (defined in either 
way) converge to that of the reflected Brownian motion on D. We remark 
that, in Theorems 5.1 and 5.6, the myopic conditioning approximation of 
reflected Brownian motion is proved for every starting point x £ D, so these 
theorems demonstrate explicitly that the symmetric reflected Brownian mo- 
tion on D is completely determined by the absorbing Brownian motion in 
D. 

We would like to point out that for the first two approximation schemes 
(i.e., random walk approximations) discussed in this paper, the proof of 
tightness of the approximating sequences requires only the weak assumption 
that D is bounded and its boundary has zero Lebesgue measure. However, 
an example given in Section 4 shows that in some domains D, reflected 
random walks do not converge to the reflected Brownian motion in D. 

The definition of the random walk on a lattice depends very much on the 
geometric structure of the state space. Myopic conditioning, on the other 
hand, can be applied to any Markov process that has a positive probability 
of not hitting a set for any fixed amount of time, if it starts outside that set. 
Hence, myopic conditioning might provide a new way of defining reflected 
Markov processes, for example, the reflected stable process that was intro- 
duced in [1], whenever the limit can be shown to exist. We plan to address 
this problem in a future work. 
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The literature on topics discussed in this paper is rather hmited. An 
invariance principle for discrete approximations to the reflected Brownian 
motion was proved in [15] in C'^-domains. Our approach in Theorem 2.4 of 
this paper is quite different from that in [15]. For results on approximations 
to the killed Brownian motion in Lipschitz domains, see, for example, [18]. 

In the rest of this introduction we give a brief review of RBM on nons- 
mooth domains, as well as present some definitions and notation, followed 
by a brief description of the approach of this paper to establish discrete 
approximations of reflected Brownian motion. 

Let n > 1 and D he a. bounded connected open set in M". Denote by m 
the Lebesgue measure in M"'. Define 

W''\D) := {/ G L\D,m):Vf G L\D,m)], 

equipped with norm ||/||i,2 := H/lb + ||V/||2, where denotes the 

U'{D,m)-\iovm. of / for p > 1. In Sections 2 and 3 we will assume that 
the boundary of D has zero volume, that is, m[dD) = 0. To prove the weak 
convergence of random walks, we need to impose the following condition on 
D: 

(1.1) C^(D) is dense in {W^''^{D), \\ ■ ||i,2). 

Here C^{D) is the space of real- valued continuous functions on D that have 
continuous first derivatives on D. Condition (1.1) is satisfied when is a 
VF"^'^-extension domain in the sense that there is a bounded linear operator T 
from (iyi'2(D), II . ||i^2) to (VF^-^^R"), || • ||i,2) so that Tf = f m-a.e. on D for 
every / G W^''^{D). This is because C^{W^), the space of smooth functions 
with compact support in M", is || • ||i^2-dense in W^''^{W^). Examples of VF^'^- 
extension domains are bounded Lipschitz domains in M", and, more gener- 
ally, local uniform domains also known as (e, (5)-domains (see [12]), defined as 
follows. We say that D is an (e, 5)-domain if 6,e > 0, and whenever x,y £ D 
and |x — y| < 6, then there exists a rectifiable arc 7 C joining x and y with 
length(7) < e~"'^|a; — y\ such that min{|a; — z\, \z — y\} < dist(z, 3-D) for 
all points z G 7. Here dist{z,dD) is the Euclidean distance between a point 
z and the set dD. "Uniform domains" can be defined as (e, oo)-domains. An 
example of a uniform domain is the classical von Koch snowflake domain. 
Every nontangentially accessible domain defined by Jerison and Kenig in [11] 
is a uniform domain (see (3.4) of [11]). The boundary of a uniform domain 
can be highly nonrectifiable and, in general, no regularity of its boundary 
can be inferred (besides the easy fact that the Hausdorff dimension of the 
boundary is strictly less than n). For any a G [n — 1, n), one can construct a 
uniform domain D C such that 7i^{U R dD) > for any open set U sat- 
isfying C/n dD 7^ 0. Here denotes the a-dimensional Hausdorff measure 
in M". 
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Under condition (1.1), the Dirichlet form {£, W^''^{D)) is regular and, there- 
fore, there is a symmetric diffusion X taking values in D associated with 
it, called the reflected Brownian motion on D. See [9] for definitions and 
properties of Dirichlet spaces, including the notions of quasi-every where, 
quasi-continuous, and so on. When D is C^-smooth, X admits the following 
Skorokhod decomposition (cf. [5]): 



where B is the standard Brownian motion in , n is the unit inward normal 
vector field on dD and L is a continuous nondecreasing additive functional 
of X (called the boundary local time of X) that increases only when X^ is 
on the boundary of D. 

Constructing a reflected Brownian motion on a nonsmooth bounded do- 
main D is a delicate problem. Fukushima [8] used the Martin-Kuramochi 
compactification D* of D to construct a continuous symmetric diffusion pro- 
cess X* on D* whose Dirichlet form is {£,W^'^{D)). The process X* could 
be called reflected Brownian motion in D, but it lives on an abstract space 
D* that contains D as a dense open set. It was proposed in [5] to refer the 
quasi-continuous projection X of X* from D* into the Euclidean closure D 
as reflected Brownian motion in D. The projection process X is a continuous 
Markov process on D and it spends zero Lebesgue amount of time on dD. 
Moreover, it behaves like Brownian motion when it is inside D. For every 
point X £ D, one can construct reflected Brownian motion in D defined as 
above so that it starts from x. The distribution of X is uniquely determined 
by the fact that its associated Dirichlet form is {£,W^''^{D)). In general, X 
is not a strong Markov process on D (e.g., this is the case when D is the 
unit disk with a slit removed). However, when condition (1.1) is satisfied, 
X is the usual reflected Brownian motion in D obtained as the Hunt pro- 
cess associated with the regular Dirichlet form {£,W^''^{D)) on D. See [2], 
Section 3 and [5] , Section 1 for more information on the history of reflected 
Brownian motion on nonsmooth domains. 

For any metric state space S and positive T > 0, let D([0, T],5) denote the 
space of all functions on [0,r] taking values in S that are right continuous 
and have left limits. The space D([0,oo),5) is a separable metric space 
under the Skorokhod topology. We refer the reader to [7] for the definition 
and properties of the Skorokhod topology and space D([0,r],5). The space 
of continuous functions from [0,T] to S, equipped with the local uniform 
topology, will be denoted by C{[0,T],S). We will use similar notation for 





t > 



6 



K. BURDZY AND Z.-Q. CHEN 



spaces of functions defined on finite time intervals. For D C M", we will use 
C^{D) to denote the space of smooth functions with compact support in 
D. 

For technical convenience, we will often consider stochastic processes 
whose initial distribution is a finite measure, not necessarily normalized to 
have total mass 1, for example, the Lebesgue measure on a bounded set D. 
Translating our results to the usual probabilistic setting is straightforward 
and so it is left to the reader. 

We close this section with a brief description of the main ideas of our 
proofs. Denote by {PP,t > 0} the transition semigroup of Brownian mo- 
tion killed upon leaving domain D. In the first two approximation schemes 
(i.e., random walk approximations), define mfc(x) = "^[^^ 2^^'" for x € D^, 
where Vk{x) is the degree of the vertex x in D^; and in the myopic condi- 
tioning scheme, define mk{dx) := l£){x)P^kl{x)m{dx) . Let X'^ be one of 
the discrete approximating processes mentioned above. Then is the re- 
versible measure for in a suitable sense. We will show that the law of 
{X^,F'^^,k> 1} is tight in the space C([0, 1],]R") or L»([0, 1], M"), and that 
any of its weak subsequential limits (Z, P) is a time-homogeneous Markov 
process that is time-reversible with respect to the Lebesgue measure m in 
D. We then show that the process Z killed upon leaving domain D is a 
killed Brownian motion in D and establish that the Dirichlet form ,J^) 
of Z in L'^[D,m) has the property that 

and 

£^{f,f)<£{f,f):=U |V/(x)|2dx iorf€W'^\D). 

JD 

We thus conclude from the following Theorem 1.1 that {£^,J=') = {£, W^''^{D)) 
and, therefore, {Z,P) is the stationary reflected Brownian motion on D. 
Discrete approximations with nonstationary initial distributions are then 
handled via those with stationary distributions. 

Theorem 1.1. LetD be a bounded domain inM" and mo be the Lebesgue 
measure in D that is extended to D by taking m,o{D\D) =0. Suppose that 
Z is a D-valued right continuous time-homogeneous Markov process hav- 
ing left-limits with initial distribution mo and is symmetric with respect to 
measure mo- Let {£^ ,J^) be the Dirichlet form of Z in L'^{D,mD)- If the 
subprocess of Z killed upon leaving domain D is a killed Brownian motion 
in D, then 



TC1W^^\D) and £^ (f, f) > £{f, f) forf^T. 
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The above theorem is essentially due to Silverstein [13], Theorems 15.2 and 
20.1. In fact, these theorems were proved in a more general context. Observe 
that the Dirichlet form for the killed Brownian motion in D is Wq''^{D)) 
and {£,W^''^{D)) is its active reflected Dirichlet space. Here Wq'\D) is the 
completion of C^{D) under the norm || • ||i,2- An accessible proof of Theorem 
1.1 can be found in Takeda [17], Theorem 3.3. 

2. Invariance principle for discrete reflected random walk. Let D be a 

bounded connected open set in M" with m{dD) = 0. Without loss of gener- 
ality, assume that the origin £ D. Let 2~'^Z" be the union of all closed line 
segments joining nearest neighbors in 2~^1^"', let be the connected com- 
ponent of 2-'=Z'^ n D that contains the point 0, and let Dk = Dln 2-'^'Z". 
For X € Dk, we use Vk{x) to denote the degree of the vertex x in D^. Let 
{X^2-2k,j = 0, 1, . . .} be the simple random walk on that jumps every 
2-2fc ^j^j^ q£ time. By definition, the random walk {X^ 2fe j J = 0, 1, . . .} 
jumps to one of its nearest neighbors with equal probabilities. This dis- 
crete time Markov chain is symmetric with respect to measure m^, where 
rukix) = H|[£l2-'=" for x e Dk- Clearly, nik converge weakly to m on D. 
We now extend the time-parameter of {Xj2-2k,j = 0, 1, . . .} to all nonneg- 
ative reals using linear interpolation over the intervals ((j - 1)2-2^=, j2-2'^) 

for j = 1,2, We thus obtain a process X'^ = {X^,t > 0}. Its law with 

Xq = X will be denoted by P^. 

For x,y G Df^, let x y mean that x and y are at the distance . Let 
Qk{x,dy) denote the one-step transition probability for the discrete time 
Markov chain {X^^_2k,j = 0, 1, . . .}; that is, for / > on and x G D^, 

Qkf{x):= f f{y)Qk{x,dy):=-^ J2 /(?/)• 

For feC^(D), define 

A/(x):= / ifiy)-f{x))Qk{x,dy) 
Jd 



^^(^^ y&Dk-.y^x 

Then - Eto CkfiX^^.2k),G%-k,j = 0, 1, . . .} is a martingale for 

every / £ C^(D), where := a{X^ ,s <t). 

To study the weak limit of {X'', > 1}, we introduce an auxiliary process 
defined by YJ^ := -'^[22'=t]2-2'= ' "^^ere [a] denotes the largest integer that 
is less than or equal to a. Note that Y^ is a time-inhomogeneous Markov 
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process. For every fixed t > 0, its transition probability operator is symmetric 
with respect to the measure on D^. Let J-j^ := a{Y^,s < t}. By abuse of 
notation, the law of starting from x G Dj. will also be denoted by P^. 

Lemma 2.1. Let D be a bounded domain in with m{dD) = 0. Then 
the laws {P'^^,k> 1} of{X'',k> 1} are tight in the space C{[0,T],W^) for 
every T > 0. 

Proof. For each fixed A; > 1, we may assume, without loss of generality, 
that O is the canonical space D([0,oo),M") and Y^^ is the coordinate map 
on Q. Given t>0 and a path uj gQ, the time reversal operator rt is defined 
by 

'uj{{t-s)-), ifO<s<i, 
uj{0), ifs>t. 

Here for r > 0, w(r— ) := lims-f^ a;(s) is the left limit at r, and we use the 
convention that a;(0— ) :=w(0). We note that 



(2.1) rtiij){s) :-- 



limrt{uj){s) = uj{t—) = rt{uj){0) and 
sJ-O 

(2.2) 

hm rj(w)(s) =uj{0) = rt{uj){t). 

Observe that for every integer T > 1, P^^. restricted to the time interval 
[0, T) is invariant under the time-reversal operator r^. Note that 

M^'^:=f{Y,')-f{Y,')- ^kf{Yt,-..) 

i=0 

is an {J-'t,t > 0}-martingale for every / G C'^{D). We have 

(2.3) f{Yt) - f{Y^) = iMf - i(M^'/ - Mfy_,)_) o rr for t G [0,r). 

For fi{x) = Xi, let M^^' := M^^f- and := (M'^'S . . . , M'^'"). With this 
notation, we have 

Y^'^ - = - \{m!}- - M(V,)J o rr for t G [0, T). 

Let 5"^ := {x G Dfc : v{k) = 2n}. For / G C^{D), by the Taylor expansion, 
for X G Dk, 



1 " d^f{x' 



+ 2 51 -Q^\riyi-^i)iyj-^j) + 0{'\)\y-x\'\Qk{x,dy). 
i,i=l * / 
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We see that, for x € Sk, 

i^d^fix 



and so 

(2.4) Ckfix) = i-A/(x)2-2fe ^ o{l)2-'"', 

while for x £ Df^ \ , 

|/:fc/(x)|<2-2^o(i). 

For 6>0, define := {x G M" : dist(x, dD) < 26}. Note that D^'"^ con- 
tains the i/n2~^-neighborhood of \ Sk and so mk{Dk \ Sk) < m{D'^ ^' 
which goes to as A; ^ oo because m[dD) = 0. Thus, for every t > 0, 

■ (/-l)2-2fe ~ ^kXi{i (/_i)2 



1=1 



[22*=t] 

{1-1)2- 



«=1 



Since is the invariant measure for the Markov chain {Y^^_2k,j = 0, 1, . . .}, 



we have 

< lim tm(D' 

n 



(!-l)2-2fe 



fc— ♦CO 



[Af'='\M^'-']t - < hm tm{D'^ 'v^)O(l) = 0. 



By [7], Theorem 7.4.1, the laws of M'^ converge weakly to that of a Brownian 
motion in M" in the space D([0,T],]R"'). Since P^j. is invariant under the 
time-reversal operator rj- when restricted to the time interval [0,T), we 
have by [10], Proposition VI. 3. 26, that the laws of {Y^,k > 1} are tight in 
D([0,r),]R"') and any of subsequential limits of {Y^,k > 1} is the law of a 
continuous process. Now by [7], Proposition 3.10.4 and [7], Problem 3.25(d) 
on page 153, we conclude that the laws of {X^,k > 1} are tight in the space 
C([0,T),M"). □ 

Let {X,P) be any of subsequential limits of (X",Pj]^J, say, along a sub- 
sequence {X"j,Pj^^_, ,j > 1}. Let To := inf{t > : Xt ^ D} and mo '■= m\D- 
Clearly, Xq has distribution m£). 
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Lemma 2.2. In the above setting, for every f G C^{D), the process 
m/ := f{Xt) — /{Xq) — 2^ /q Af{Xs) ds is a -square integrable martingale. 
This, in particular, implies that {Xt,t < rD,P} is a Brownian motion killed 
upon leaving D, with initial distribution mo and infinitesimal generator 
f A. 

Proof. The lemma follows easily from the invariance principle for ran- 
dom walks, but we sketch an argument for the sake of completeness. Recall 
the definition of {X^,Y'^,x € D) and {Qt,t > 0}, the cj-field generated by 
X''. For / G C^iD), there is a ko such that 

supp[/] C G i:):dist(2;,D") > 2-2^0 

Thus, by (2.4), 2^^£^./ converges uniformly to ^^A/ on D. 

Without loss of generality, we take the sample space of X^ and X to 
be the canonical space C([0, By the same argument as that in the 
last paragraph on page 271 of [16], we can deduce that {m/ ,t > 0} is a 
P-martingale. □ 

Lemma 2.3. Let D be a bounded domain in and fix k> 1. Then for 
every j>\ and f £ L'^{D,mi^), 

if - Q^k f, f)LHD,mk) < J(/ - Qkf^ f)L^{D,mk) ^ 2i(/ - Qkf, f)L\D,ni^)- 

Proof. Note that Qfc is a symmetric operator in L'^{D,mk)- So for 
feL\D,mk), 

{Qkf - Qkf^ f)v^(D,m^) = {f - Qkf^ Qkf)L'2(D,mt,) 

< (/ - Qkfif)L^iD,mk)- 

Moreover, 

{Qk9,a)L2{D,mk) = {Qkg,Qkg)L^{D,mk) ^0- 

We apply the last remark to g := f — Q\f to obtain 

{Qkf - Qkfi f)L'2{D,mk) = {Qkf - Qkf: Qkf)L^{D,mk) 

= {Qkf - Qkf-: f)L^(D,mk) 

-{Ql{f-Qlf),f-Qlf)LHD,m,) 
< (Qkf - Qkf-: f)L2{D,mk)- 

Suppose the following holds for some j > 2: 

{Qff - Q^^^^f, f)L'^{D,mk) - (QT f - Qff^ f)L^{D,mk) 

for every i < j. 



DISCRETE APPROXIMATIONS TO REFLECTED BROWNIAN MOTION 11 



Then 

Wk J~^k J^j)L^iD,mk) 

= {Ql^'''\Qlf)-Ql'{Qlf),Qlf)LHD,m,) 

< iQl^'-'\Qlf)-Qf'-'\Qlf),Qlf)LHD,m,) 

= {Qk f ~ Qk'''^^^ f^f)L■2{D,mk)■ 
This proves by induction that, for every i > 1, 

{Q'k f - Qk^^""^ f ^ f) L'^iD^nik) - ('Sfc^* f ~ f ^ f) L-^{D,mk)- 

It follows that 

(Qi7 - f)LHD,m,) < (/ - Qlf^ f)LHD,m,) for every i > 1, 

and so 



(2.5) ^ 

Since Qk is a symmetric operator in L^{D,'mk), we have 

(/ - Qkf, f)L^{D,mk) = if - Qkf, f)L2{D,mk) + (/ " Qkf, Qkf)L^{D,rrn,) 
< 2(/ - Qkf,f)L^(D,rrn,)- 

This and (2.5) prove the lemma. □ 

We will say that "Z^ is a Brownian motion running at speed 1/n" if Znt 
is the standard Brownian motion and we will apply the same phrase to the 
other related process. 

Theorem 2.4. Let D be a bounded domain in M" with m{dD) = and 
satisfying the condition (1.1). Then for every r>0, the laws o/{X^,P^J 
converge weakly in C([0, TjjM") to a stationary reflected Brownian motion 
on D running at speed 1/n whose initial distribution is the Lebesgue measure 
in D. 

Proof. Fix T > 0. Let P) be any of the subsequential limits of 
(X*^,P^J in C([0,r],M"), say, along (X'^^P^,,). Clearly, X is a time- 
homogeneous Markov process with transition semigroup {Pt,t > 0} that 
is symmetric in L'^{D,dx). Let {Pt,t G 2~'''Z_|_} be defined by Pj'f{x) := 



12 K. BURDZY AND Z.-Q. CHEN 

^xifi^t)]- For dyadic t > 0, say, t = 30/2^^° and / G C^{D), we have, by 
Lemma 2.3 and the mean-value theorem, 

-(/ - Ptf,f)L'2(D,dx) 

1 fc, 

= ^ ^.1™ (/ - V, f)L^D,mk^ ) 

= — Um(/-Qif fJ)LHD,m,.) 
92/00 

< linisup— io22'=^-2'=«(/ _ Q,JJ)^,^^^^^ ^ 

J— >oo Jo 



hmsup2(2-")'^^i- ^ 5] (/(^)2_;(^)/(y)) 
hmsup2(2-«)fc^-^ ^ ^ (/(^)2_/(^)/(y)) 



+ hmsup2(2-)'=.J- ^ ^ 

ySDfc^, xeDfe^. : y^x 

:hmsup2(2-)'=.i- ^ (/(^)_/(y))2 

x,yf^Dk,:y'--*x 



-f 

2nJD 



|V/(x)|2dx. 



Let {£,J') be the Dirichlet form of X, or equivalently, of semigroup {Pt,t> 
0}. That is, 

= (/ G L2(D, dx) : sup ^(/ - PJ, f)mD,dx) 
I i>0 I 

= l™^(/--Pt/,/)L2(B,d^.) <CX)|, 

£{f,f)=Snp];{f-PtfJ)L2(D,dx) 
t>0 t 

= l^l\if - Ptf^ f)L2(D,dx) for / G J^. 



Then for feC^iD), 



£{f, /) = lini i(/ - Ptf, f)L2^n,dx) ^ ^ I W( 



x)\'^ dx. 
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This shows that / G V{£). As C^{D) is dense in {W^'^{D), \\ ■ ||i,2), we have 
c:r and 

£{f,f)<^ f \Vf{x)\^dx for every feW'^\D). 
2n Jd 

This, Lemma 2.2 and Theorem 1.1 imply that J- = W^''^{D) and 

£{fJ) = ^l |V/(x)|2(ix for /EW^i'2(D). 
2n Jd 

We deduce then that X is a stationary reflected Brownian motion on D 
running at speed 1/n. This proves that converge weakly on C([0, T], M") 
to the stationary reflected Brownian motion on D. □ 

Remark 2.5. By [7], Proposition 3.10.4, under the assumptions of The- 
orem 2.4, the stationary laws Pj^^ of the step processes defined at the 
beginning of this section converge weakly in D([0,T],]R'^) to the stationary 
reflected Brownian motion on D running at speed 1/n, for every T > 0. 

3. Continuous-time reflected random walk. Let D be a bounded domain 
in with m{dD) = and let be defined as in the previous section. But 
in this section, will be the continuous time simple random walk on Dj^, 
making jumps at the rate 2~'^^ . By definition, X^ jumps to one of its nearest 
neighbors with equal probabilities. This process is symmetric with respect 
to measure m^, where mk{x) = "^[^^ 2~^" for x G D^. Note that converge 
weakly to the Lebesgue measure m on D, and recall that, for x,y £ D/^, we 
write X '^y ii x and y are at the distance . The Dirichlet form of X^ is 
given by 

^'(/'/) = ^ E 2-(«-2)^(/(x)-/(y))l 

Lemma 3.1. Let D be a bounded domain in R" with m{dD) = 0. Then 
forfeC\D), 

(3.1) lim g'^{f,f) = l- I \Vf{xtm{dx). 

Proof. For 5 > 0, define := {x G M" : dist(x, dD) < 6}. As m{dD) = 
0, for every e > 0, there is > such that m{D^) < e. Take integer kQ>l 
large enough so that 2~^^^ < 5. Recall that for x G Dj^, Vk{x) denotes the 
degree of vertex x in the graph D^. Define '■= {x G : v{k) = 2n}. Then 
for k > ko, {D \ D^) n 2^^!/' C Sk- As 

^'(/' /) = E ( E 2"("-2)'=(/(x) - f{y)f\ , 
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we have by the mean-value theorem that 

<0{m{D^)). 



hm sup 

fc— >oo 



2n Jd\d^ 

Taking 5 yields the claim that \\Ta.k_ao£^{f, f) = Id | V/(x)pm((i2;) 
for f£C^(D). □ 

Let P^^ denote the distribution of {X^,t > 0} with the initial distribution 
ruk. 

Lemma 3.2. Assume that D is a bounded domain in M" with m{dD) = 
0. For every T > 0, the laws of stationary random walks {X^,P!^^,k > 1} 
are tight in the space D([0,r],L)) equipped with the Skorokhod topology. 

Proof. For constant T > 0, let r^ be the time-reversal operator from 
time T for X'' [see (2.1) for its definition]. Note that P^^-a.s., X'' is con- 
tinuous at time T and P^^. is invariant under rx- For / G C^(D), we have 
by [6], Lemma 3.5 and (3.6) the following forward-backward martingale de- 
composition for f{Xf ). For every T > there exists a martingale M'''^ such 
that 

(3.2) fix,') - f{Xl) = iMf - i(M^i - Mf/^^_,) J o r^+i 

for t G [0,T]. 

Note that the symmetric jump process X' has Levy system {N^ {x , dy) ,t) , 
where, for x G D^^ 

N\x,dy) = ^ E 2-^--^)H-H^y}{dy) = ^ ^ ^^H^y^idy). 
yeDk ■■ y^x y&Dk ■■ y^x 



Thus, 

{M'^f)t = [\fiX^) - f{y)fN\xldy)ds 
Jo 

= ^f E 2''^{f{X'!)-f{y)fds. 
y&Dk-.y^Xk 

Recall that v{k) is the degree of vertex x in the graph D^- Taking /(x) = Xj, 
i = 1, . . . , n, we have, for every k>l and t > s > 0, with M*^'* = M^'^^ , 

Y.m'''')t - {M''')s) < — / Vk{X^)2^'^2~^'^dr <t-s. 
i=i ^'^ -^^ 

This implies that the sequence {T,i=i{M'''')t,k> 1} is C-tight in D(M) in 
the sense of [10], Proposition VI. 3. 26. Hence, by [10], Theorem VI. 4. 13, the 
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laws of {(M'''\...,M'''"),fc > 1} are tight in D([0,r],M"). As converges 
weakly to m, and is invariant under vt+i for every k> 1, we conclude by 
[10], Theorem VI.3.21, that the laws of{X'',k>l} are tight in D([0, T], M") 
and, hence, on D([0, oo), D). □ 

Theorem 3.3. Let D be a bounded domain in M" with m{dD) = and 
satisfying the condition (1.1). Then for every T > 0, the stationary random 
walks on converge weakly in the space D([0, T],Z)), as k^oo, to 
the stationary reflected Brownian motion on D running at speed l/n, whose 
initial distribution is the Lebesgue measure in D. 

Proof. Let {Z,P) be any of the subsequential limits of (X'^,P^^), say, 
along X'^J . Clearly, Z is a time-homogeneous Markov process under P and its 
transition semigroup {Pt,t > 0} is symmetric with respect to the Lebesgue 
measure m on D. By a similar argument as that in the proof of Lemma 
2.2, the process Z killed upon leaving D is a killed Brownian motion in D 
with speed l/n. Let {S,J-) be the Dirichlet form associated with Z, and let 
{Pf ,t > 0} be the transition semigroup for X'^. As X^^ converge weakly to 
Z, we have, for every / G C'^{D) and t > 0, 



hm i(/-Pfv,/)= lim Ie 

J— »00 t J^OO t t'm,, 

= iEp[/(Zo)(/(Zo) - f{Zt))] = -^{f- Ptf, f). 
Thus, for / G C^p), by Lemma 3.1, 



By assumption (1.1), C^{D) is dense in the Sobolev space VF^'^(D) with 
respect to norm || • ||i^2- It follows that T D W^'^{D) and 



£if,f) 



snphf-PtfJ) 




< 
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Define 

= ^ f V/(x) • Vg{x)m{dx) for f,ge W''\D). 

Note that Ty^'^(Z))) is the Dirichlet form for the reflected Brownian 
motion on D running at speed 1/n. On the other hand, as we have ob- 
served at the beginning of this proof, the process Z killed upon leaving D 
is a killed Brownian motion in D with speed 1/n. Therefore, according to 
Theorem 1.1, {S,T) = {8^ ,W^''^{D)). In other words, we have shown that 
every subsequential limit of is reflected Brownian motion on D with ini- 
tial distribution being the Lebesgue measure on D and with speed 1/n. This 
shows that converges weakly on the space D([0,oo),-D) to the stationary 
reflected Brownian motion X on D running at speed 1/n. □ 

4. Examples. All the results in the previous two sections apply to any 
bounded domain D that satisfies the condition (1.1) and whose boundary 
has zero Lebesgue measure. As we noted in Section 1, bounded uniform do- 
mains have such properties. Bounded Lipschitz domains and bounded non- 
tangentially accessible domains are uniform domains. Although the Haus- 
dorff dimension of the Euclidean boundary of any uniform domain D C is 
strictly less than n and, thus, dD has zero Lebesgue measure in M", dD can 
be highly nonrectifiable. The classical "von Koch snowflake" planar domain 
defined below is such an example. 

To define the von Koch snowflake, start with an equilateral triangle Ti. 
Let / be any of its sides. Add an equilateral triangle such that one of its 
sides is the middle one third of / and its interior does not intersect Ti . There 
are three such triangles; let T2 be the closure of the union of these three tri- 
angles and Ti . We proceed inductively. Suppose / is one of the line segments 
in dTj. Add an equilateral triangle such that one of its sides is the middle 
one third of I and its interior does not intersect Tj. Let Tj+i be the closure 
of the union of all such triangles and Tj. The snowflake DvK is the interior 
of the closure of the union of all triangles constructed in all inductive steps. 

It is elementary to check that DvK is a nontangentially accessible domain 
and so a uniform domain or an (e, oo)-domain (see Section 1 for definitions). 
It is also well known that the Hausdorff dimension of dD^K is wf- Hence, 
the 2-dimensional Lebesgue measure of dD^x is 0. We conclude tnat all re- 
sults stated in the previous two sections, in particular. Theorems 2.4 and 
3.3, apply to the von Koch snowflake. 

Without some domain regularity conditions, the results in the previous 
two sections do not have to be true. Here is a counter-example. Let 

U^ = {{x,y) £ (0,1)^:|2;- j2"''| < e or \y-j2~^\ <e for some j G Z}. 

We choose > so that \U^''\ < 2'^^^ and let U = Ufc>i Ul'' . Note that 
i7 is a bounded open connected set with Lebesgue area less than 1/2. Let 
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Dk be defined as in the previous section, relative to D = (0, 1)^. Note that 
for every A; > 1, Dk GU, so the sets defined relative to U are the same 
as those defined relative to D = (0,1)^. It follows from Theorems 2.4 and 
3.3 that processes X'^, with discrete and continuous time, defined relative 
to Dk, converge weakly to the reflected Brownian motion in [0,1]^. Since 
m{U) 7^m((0, 1)^), the reflected Brownian motion in U has a different dis- 
tribution than the reflected Brownian motion in (0,1)^ and it follows that 
the conclusions of Theorems 2.4 and 3.3 do not hold for U . 

We would like to emphasize the fact that the approximation scheme for 
reflected Brownian motion developed in the next section works for any 
bounded domain. 

5. Myopic conditioning. Throughout this section, D C M" is a bounded 
connected open set and X = {Xt,t > 0, P^:, x G M"} is a Brownian motion in 
M". Let = {Xj^^t >^,Y'^,xeD] he the Brownian motion in R" killed 
upon leaving the domain D. For each /c > 1, we deflne a myopic process X^ as 
follows. For t E [0,2"*^], let X^ be X^ conditioned not to leave domain D by 
time . Suppose that X^ is now defined on the time interval [0,j2^'^]. We 
deflne X^^+s * ^ (0' 2"''] to be a copy of X^ conditioned not to leave 
domain D by time 2~'^, starting from Xjg-fe, but otherwise independent of 

{X^,t£ [0, j2-*^]}. The law of X^ with X^ = x wiU be denoted as and 
the mathematical expectation under P^ will be denoted by E^. 

Theorem 5.1. For every hounded domain D in M" and for every xq S 
D, the processes X^ under P^^ converge weakly to the reflected Brownian 
motion Y on D starting from xq in the space C([0, 1],-D) as k ^ oo. 

To prove the above theorem, we introduce auxiliary processes in which 
pieces of conditioned Brownian paths are replaced with line segments. More 
precisely, we let {Y'',t G [0, 1]} be constructed from {Xj2_k,j = 0, 1, . . . ,2*^} 

by linear interpolation over the intervals {{j — l)2~^,j2~^) for j = 1, . . . , 2^^. 
Note that {Y!^^_^.,j = 0, 1, . . . , 2^^} is a Markov chain with one-step transition 
probability Q^, where, for x £ D and k>l, Qk{x,dy) is the distribution at 
time of X^ that starts from x and is conditioned not to leave D by time 
2~'^. In other words, if we let {PP,t > 0} denote the transition semigroup 
for the killed Brownian motion X^ , then for any nonnegative Borel function 



With a slight abuse of notation, the law of Y^ with Yq = x will also be 
denoted as P^ and the mathematical expectation under P^ will be denoted 



/on D, 
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by E^. Let mk{dx) := l£){x)P^k^(x) dx. Observe that 

(5.1) {Qkf,g)L^(D,mk) = if^Qk9)L2{D,mk) foT f , g > OU D , 

and so is a reversible measure for Markov chain {Y^^_f,,j = 0, 1, . . . , 2'^}. 
Let m o denote the Lebesgue measure on D that is extended to M"" by letting 
mD(M"' \ D) = 0. It is clear that converge weakly to m^) on D. 

We will show that converge weakly to reflected Brownian motion on 
D and then use this fact to establish Theorem 5.1. 

Lemma 5.2. Suppose that either (i) /i^ = for every k>l; or (ii) 
{/J'kjk > 1} is a sequence of measures on D with supi^^i iJ,k{D) < oo and 
HkiD \ K) = for some compact subset K of D and all k>l. Then the 
laws o/{y'=,P^^,/c> 1} are tight in the space C{[f),l],W). 

Proof, (i) We first prove the lemma under condition (i). For nonnega- 
tive / e C^{D), by Ito's formula, 



Qkf{x) 



-E.[/(Xf.,)] 



P^_,lix 



fix) + -E,. 



r A/(Xf)ds 

JO 



>f{x) 



IIA/II 



Fix k>l. Let C/*^ = a{Yt^_^,i < j). For nonnegative / G C^{D), let 



IIA/II. 



-. Then we see from the above that 



Elj/(y(5+,)2-.) + Af{j + l)2->'\g^,.,] = Q,f{Y^^.,) + Af{j + 1)2-^ 



In other words, {/(y^2„fe) + ^/i2-^g|2-Ji=o,l,...,2 
Pj^^-submartingale. Moreover, for every e > 0, 



fe is a nonnegative 



fc— >oo 



0-1)2^ 



< lim 2'' [ P:r(|^2-'= - ^ol > e and 2''' < td) dx 

< lim 2^= / P^{\X2-k-Xo\>e)dx 



0. 
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Thus, by [16], Theorem 1.4.11, the laws of {Y^,P'^^,k > 1} are tight in 
C([0,1],M"). 

(ii) Now assume that {fJ-k,k > 1} is a sequence of measures on D satisfying 
condition (ii). Since every compact set has a finite covering by open balls, 
we may and do assume, without loss of generality, that K C B{xQ,rQ) C D. 
Define 6o := dist{B{xQ,rQ), D'^)/2 and B := B{xo,ro + do). Recall that X is 
Brownian motion in M". By Lemma II. 1.2 in Stroock [14], there is a constant 
Co > so that 

(5.2) p(^sup|X,-Xo| >r^ <coexp(^--^^. 
This, in particular, implies that 

(5.3) P4^1(x) > 1 — coexpf for every x G -B. 

V Cot J 

Let denote the semigroup for Brownian motion killed upon exiting B. 

Without loss of generality, we take the sample space of and X to be 
the canonical space C([0, 1], M"). For a; G C([0, 1], M") and /9 > 0, we define 
the oscillation of uj over time interval [so,to] by 

oscp[so,to]{uj) ■= sup \u;{s) — u}{t)\. 

s,te[so,to] : \t-s\<p 

We also define 

To :=inf{i > 0:w(t) ^ B}. 

Sometimes we will add a superscript to make the underlying process explicit, 
for example, we may write P(osc^[so,to] > £)• 

Since all n^s are supported in a compact set and have uniformly bounded 
mass, standard theorems (see, e.g., [16], Theorem 1.3.1) show that the laws 
of {y^,P^^, A; > 1} are tight in the space C([0, 1],]R") if and only if, for every 
e>0, 

limsupP^ (oscp[0,l] >e) =0. 

PlO k>l 

So it suffices to show that, for every e > and 6 > 0, there is p > and 
N>1 such that 

(5.4) P^^ (oscp[0, l]>e)<6 for every k>N. 

Fix an arbitrarily small S > 0, and let to = jo'2^° be a dyadic rational in (0, 1) 
such that coexp{—-^^)supi^yifj.k{D) < 5/4. Define 

Cfc := y^ - Co exp — —j j 
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and note that Cfc ^ 1 as /c — > oo . We have 

P^^ (oscp[0, 1] > e) < P;;^ (ro < to) + Pj^, (oscp[0, 1] > e and tq > to)- 

Observe that B is convex. For k> ko, it follows from the definition of the 
process and (5.3) that 

P'.(^o<to) 

= EPM.(^oG((j-l)2-^i2-^•]) 



(5.5) 



2H0 



n 



^2-^^^-^?2-^^ ) 



Xg.^ G S for i = 0, . . . , j - 1 but ^ S 



< Cfc 5] P/,, (X^^.fe G 5 for i = 0, . . . , j - 1 but ^ S) 

So 



< CfeCo exp 

< Ck6/4. 



coto 



On the other hand, by the definition of Y'' and (5.3) again, conditioned on 
{tq > to}, the law of restricted to {0,to] is dominated by Ck times that 
of linear interpolation of X^ at times j2~^, with initial distribution fi^- 
Thus, for any e > 0, one can make to> smaller, if necessary, so that for all 
sufficiently large k, 



pfc 



{tq > to and oscp[0,to] > £/2) 



(5.6) 



< CfcPp,, {t^ > to and osc^ [0, to] > e/2 



< Cfc^fc(D)Po(osc; [0,to] > e/2), < 5/4, 

where Pq denotes the law of Brownian motion X starting from the origin. 

For each fixed y G B{xo,ro), let x 1-^ ■ipi:{y,x) be the density function for 
the distribution of Xf^ under Py restricted on the event that {X^_^ G 



B for j = 1,2, . . . Clearly, x ^ ^fc(y,x) is a bounded function on D 

that vanishes outside B, and as k]oo, ^fc(y,x) decrease to (j){y,x), the 
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probability density function for the killed Brownian motion in B at time 
starting from y G i?(xo,To). There is a constant ai > such that 

sup i^kiu, 3s) dx < aiPf^l{x) dx < aimk{dx) 
yeB{xo,ro) 

for every k>l such that < to- As sup^>i fJ-kiD) < oo, there is a constant 
02 > 0, independent of A: > 1, such that the distribution of X^^ under P^^ 
restricted on the event that {X^_^. S -B for j = 1, 2, . . . , is dominated 
by a2mk{dx). We obtain, 



Pu,(7-o > to and oscpfto, 1] > e/2) 



(5.7) 



>D ^ Pxf (oscp[0, 1 - to] > e/2); 



for j = l,2,...,2% 



, (24o+l)/(2'=io) pfc / rp, 11^^ /ON 

<c^ a2Pmfe(osCp[U, IJ > e/2j. 

We have already proved in part (i) of the proof that the laws of {Y^, P^j. ) ^ ^ 
1} are tight in C([0, 1],1R"), so there is TV > 1 such that 

(5.8) P^^ (oscp[0, 1] > e/2) < <5/(4a2) for every A; > iV. 

Combining (5.5)-(5.8), we obtain for large k, 

P^Josc,[0,l]>e) 

< Pp, (to < to) + Pp, (oscp[0, to] > e/2 and To > to) 

+ P^^, (oscp[to, 1] > e/2 and tq > to) 

<c.V4 + 5/4 + cf*»+^)/(^'^'*")a2<5/(4a2). 
This proves (5.4) because limfc^ooCfc = 1- □ 

Consider a sequence of finite measures fik on D, k >1, that satisfies the 
assumptions of Lemma 5.2(i) or (ii) and converges weakly to a finite measure 

By Lemma 5.2, the laws of {y'', P^^ , A; > 1} are tight on C([0, 1], M"). Let 
(Y, P) be any of its subsequential limit, say, along a subsequence {y^^ , Pj]^ 
1}. Clearly yo has distribution fi. 

Lemma 5.3. In the above setting, for every f G C^{D), m/ := f{Yt) - 
/(Yo) — ^ Jq Af{Ys) ds is a -square integrahle martingale. This in partic- 
ular implies that {l(,t < t/),P}, with td := inf{t > 0:1^ ^ D}, is the killed 
Brownian motion in D with initial distribution ^. 
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Proof. Recall the definition of (y^,Pj;,x G D) and denote by {Gt,t> 
0} the cr-field generated by Y''. For / G C^{D), define 

Ckf{x)= I {f{y)-f{x))Qk{x,dy). 

JD 

Then {/(y.^O - ELi Ckf{Y^2.,),G%-,,3 = 0, 1, . . . , 2^ P^!} is a martingale 
for every / G C^{D). For / G C^{D), using the Taylor expansion, we have 

= 2'^ J^(vf{x){y-x) 

+ ^ E f/i^^. iVi - Xi){yj -Xj) + 0{l)\y - xAQk{x,dy) 

i,j=i * J 

+ 0(1)E,[|X2-. -Xo|3;2-^'<rB]j. 

This converges uniformly to ^A/(x) on as k ^ oo. 

Without loss of generality, we take the sample space of Y'' and Y to be 
the canonical space C([0, 1],]R"). Then by the same argument as that in the 
last paragraph on page 271 of [16], we can deduce that {M^ ,t > 0} is a 
P-martingale. □ 

Lemma 5.4. Let D be a bounded domain in M" and fix k> 1. Then for 
every j >l and f G L?'{D,mk), 

if - Qif, f)L^{D,mk) <j{f- Qkf, f)L^{D,mk)- 

Proof. The idea of the proof is similar to that for Lemma 2.3. Note 
that Qk is a symmetric operator in L^{D,mk)- So for / G L'^{D,mf^), 

{Qkf - Qkf, f)L'2{D,mk) = if - Qkf, Qkf)L^{D,mk) ^ if - Qkf, f)L^{D,mk)- 



We have 



{Qkg,g)L^{D,7nk) = iP2-''9,9)L2{D,dx) = J^P^-k-igixf dx > 0. 
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We apply the above observation to g := f — Qkf to obtain 

{Qkf - QkfJ)L^iD,mk) 

= {Qkf - Qkf, Qkf)L-2(D,mk) 

= {Qkf - Qkf^ /)L2(D,mfc) - {Qk{f - Qkf),f - Qkf)L'2{D,mk) 
< {Qkf - Qkf, f)L2(D,mk) ■ 

Suppose the following holds for j >2: 

{Qif - Q^k^f, f)LHD,m^) < (<3i"V - Qif, f)LHD,m^) for every i < j. 
Then 

{Qi^'f - Qi~''fJ)LHD,m,) = {Qi'\Qkf) - Qi{Qkf),Qkf)LHD,m,) 

< {Ql\Qkf)-Ql\Qkf),Qkf)L^^n,m,) 
= {Qif-Qi^'f,f)LHD,^.y 

This proves by induction that, for every i > 1, 

{Qkf - QV'^f,f)L'2{D,mk) < {QV^f - Qkf,f)L'2(D,mk)- 

It follows that 

{Qkf - Q'k^f^ f)L^D,mk) <{f - Qkf, f)L^{D,mu) fo^ ^very i > 1, 
and so 

3 

(/ - Qkf^ f)L'^{D,mk) = ^{QV^f ~ Qkf, f)L'^{D,mk) 

i=l 

< j{f - Qkf, f)L'^{D,mk)^ 

which proves the lemma. □ 

Theorem 5.5. Let D he a hounded domain in M". The processes 
under converge weakly to the stationary reflected Brownian motion Y 
on D in the space C([0, 1],M'") as k ^ oo, where Yq is distributed according 
to the Lehesgue measure on D. 

Proof. Let {Y,P) be any of the subsequential limits of (F'^^P^^) in 

C([0,1],M"), say, along {Y^^,P2,.). Clearly, Y is a time-homogeneous Markov 

process with transition semigroup {-Pj , t > 0} that is symmetric in L^{D,mD) = 
L^{D,dx). Let {Pt^,te 2-^'Z+} be defined by Pt^f{x) := E^[/(yj'=)]. As 
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,Pmkj) converges weakly to {Y,P) in C([0, 1],]R"), we have for dyadic 
rational t > 0, say, t = jo/2^° and (/> G C{D), 

lim {Pt^ct>A)LHD,m,.) = lim eJ;^ ,[</.(yo'^0</'(>^^)] 

(5.9) 

= E[</.(yo)0(i^t)] = m(/')L2(D,,n)- 

Recall from (5.1) that m/j is a reversible measure for operator Pj^ for f G 
2~'^Z+ when k is sufficiently large. For t G 2~'^Z_|_, the operator Pj^ has 
density function p^{t,x,y) with respect to the reversible measure m^. The 
density function p^{t,x,y) is symmetric in {x,y). Note that P/'T = 1 and so 

/ p'^{t,x,y)mk{dy) = 1 for every x G D. 
JD 

For a dyadic rational t £ {0,1], g £ Lp'{D,'mk) and large k, 
{g-Pt^g,g)L2 

{D,mk) 



{g{x) - Ptg{x))g{x)mk{dx) 

D 



(5.10) 



D\^^^ ~ Jd^ (*,a:^>y)f(y)"ifc(c^y)j5'(2;)mfc(d2;) 
{9{xf - 9{x)g{y))p^{t,x,y)mk{dx)mk{dy) 

DxD 



h / {g{xf - g{x)g{y))p^{t,x,y)mk{dx)mk{dy) 

'DxD 



+ h {g{yf - g{x)g{y))p'^{t,x,y)mk{dx)mk{dy) 
JdxD 



1 / „^„,^^2„fc 
2 



ig{x) - giy)) P'it,x,y)mk{dx)mkidy) 

DxD 



lBijg{Y,')-g{Y,')f]>0. 



For later reference, we record the following consequence of the last formula, 

(g - Qkg,g)L2(D,m^:) 
(5.11) = {g-P^^,g,g)LHD,m,) = '^Bi^[{g{Yo') - giY,'^,)f] 

= BMXo) - ff(^2-0)'; 2"'^ < td], 
where X is Brownian motion in and td ■= inf{t > 0:Xt ^ D}. By (5.10), 
we have for dyadic t > and large A; > 1 , 

< (5 - P2tg,g)L^D,mk) = (.g^g)L^D,mk) - {P2tgig)L2(D,mk) 

= ig,g)L2{D,mt:) ~ iPtg^Ptg)L2{D,m^)i 
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and thus, 



I P^nW'^ < llnl|2 < ||„||2 

y\\L'2{D 



Every / G L'^{D,m) can be approximated in L^(D,m) by continuous func- 
tions on D, so we have from above and (5.9) that 

^^{Pt' f,f)L^{D,mk^) = {Ptf,f)L^(D,m) ^r every / G L^{D,m). 
(5.12) ' " 

Let Z be the reflected Brownian motion on D, obtained from quasi-continuous 
projection on D of the reflected Brownian motion Z* on the Martin-Kuramochi 
compactification D* of D. Recall that Z behaves like Brownian motion in 
D before Z hits the boundary dD. For / G VF^'^(Z)), it admits a quasi- 
continuous version on D* and, therefore, f{Zt) = f{Z*) is well defined. 
Moreover, f{Zt) has the following Fukushima's decomposition (see [8] or 
[9]): 

f{Zt) - f{Zo) = Mi + nI for every t > 0, 

where is a continuous martingale additive functional of Z* with quadratic 
variation process {M^)t = /q | V/(Zs)p (is and is a continuous additive 
functional of Z* having zero energy. In particular, we have 

limiE„[(/(Zi)-/(Zo))2] = limiE„[(M/)2]= / \Vf{x)\''dx. 

These observations, (5.11)-(5.12) and Lemma 5.4 imply that for / G W^''^{D) 
and dyadic t > 0, 

1 Ik 

-{f - Ptf, f)mD,dx) = 7 if -Pt' /, f)mD,mk.) 

= — lim (/-Qif^ /,/) 

< limsup — io2'=^-'=« (/ - Qfc^ /, f)LHD,m, . ) 
j— >oo Jo ^ 



= hmsup2'=^-iE^[(/(Z2-., ) - f{Zo)y; 2~''^ < td] 
<limsup2>'^-'Bm[ifiZ^~,,)-fiZo)f] 

= lJ^\Vf{x)\'dx. 
Let {£,J^) be the Dirichlet form of Y in L'^(D,mD) =L'^{D,dx). Then for 



£{f, f) = lim i(/ - Ptf, f)mDAx) <\j^ I V/( 



x)]"^ dx. 
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This shows that / € JF. So we have W^'^{D) C J- and 

Sif,f)<U |V/(x)|2(ix for every /eVFi'2(D). 
Jd 

If follows from Lemma 5.3 that {Yt,t < t^} is a Brownian motion in 
D with initial distribution m£). Therefore, we have by Theorem 1.1 that 
J^ = W^'^{D) and 

£{f,f) = U' |V/(x)|2dx iorfGW''\D). 
Jd 

We deduce that y is a stationary reflected Brownian motion on D. This 
proves that Y'^ converges weakly on C([0, Ijjffi*^) to the stationary reflected 
Brownian motion on D. □ 

Theorem 5.6. Let D be a bounded domain D inM". Let {fik, k>l} be a 
sequence of measures on D that converges weakly to fj, with sup^>]^ fJ-kiD) < oo 
and that there is a compact set K C D such that Hk{D \ K) = for every 
k>l. Then the processes Y'^ under 'P^^, converge weakly to the reflected 
Brownian motion Y on D with initial distribution fi in the space C([0, IJjM") 
as k ^ oo. 

Proof. Using a partition of K with a finite covering of open balls if nec- 
essary, we may assume that K C B{xq, tq) C D. Let 5q = dist(i?(xo, ro),D'^)/2 
and define B : = B{xo, tq + 5o)- So the distance between ball B and D'^ is 6o. 
Define 

TB:=mi{t>0:Yt^B} and :=m^{t>0:Y^'' B}. 

According to Lemma 5.2, the laws of {{Y^ ,P^^),k > 1} are tight in the 
space C([0, 1],M"). Let {Y,P) be any of the subsequential limits of {Y^,P^J, 

say, along {Y i , P^^^. , ) . It suffices to show that the finite-dimensional distribu- 
tions of (y, P) are the same as those of reflected Brownian motion with initial 
distribution ^. For this purpose, take Q <ti <t2 < ■ ■ ■ <tN and nonnegative 
fi e Cfe(M"') for i = l,2,...,iV. We know from Lemma 5.3 that, under P, 
{Yt,t <td} is the killed Brownian motion in D with initial distribution ^i. 
Thus, for every e > 0, there is a dyadic rational to = 3o'^~^° £ (0,ti/2) such 
that 

P{tb < to) < e. 

h k ■ 

Since [Y^ ,Pf^^f.,) converge to (y, P) weakly, we have 
(5.13) limsupPjt < to) < P{tb < to) < e. 
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By (5.3), there is a constant cq > such that 



P<^l(x)>l-coexp(-^ 

Cot 



for every x £ B. 



Let PP denote the semigroup for Brownian motion kihed upon exiting D 
and be the Brownian motion X killed upon leaving ball B. The ball B 
is convex so for any bounded nonnegative function / on D and k > ko, 

= E^^ [/(y^o); >;'2-^ G S for J = 1, 2, . . . , 2H0] 



(5.14) 



E 



n 



1 



P,^_,,1{X .„,.) J P„f.l{X 



Xf^^,£B (or j = l,2,...,2%-l 
< CfcE^, [/(X,^); G S for J = 1, 2, . . . , 2%] 



where 



Cfc := 1 -coexp 



2% 
' 2co 



as /c ^ 00. 



For each fixed y G B{xQ,rQ), let x 1— > ipkiy,^) be the density function for 
the distribution of X^ under Py restricted on the event that {X^_^., G 

B for j = 1, 2, . . . , 2'^to}- Clearly, x ^ ipk{y, is a bounded function on D 
that vanishes outside B, and as A; t 00, '4)k{y,x) decrease to {to,y,x), the 
probability density function for the killed Brownian motion in B at time Iq 
starting from y. We clearly have 



(5.15) 
Let 



Et[fiYh-,4>to]>E,^[f{Xf)]. 



(t)k{x):= I p^{to,x,y)fik{dy) and (j){x) := I {to,x,y)n{dy) 



D 



D 



and denote {(j)m){dx) := cj){x)m{dx) and {(j)jm){dx) := cj)j{x)m{dx). By (5.13)- 
(5.15), we have, for each j > jo. 



4>m 



N 



to) 



.i=l 



N 

eX{\\h 

i=l 



<e1 



N 



.i=l 



— ^^j ^(bi- ■ m 



N 



i=l 



N 



i=l 
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Let Z be the reflected Brownian motion on D with initial distribution [i. 
Since <^ and are bounded continuous functions with compact support 
it fohows from Theorem 5.5 that 

N 



and 



hm E^' 



Taking j ^ oo, we have 



.1=1 



N 



to) 



to' 



N 

Y[ fii^ti-to) 
U=l 



.1=1 



N 



to J 



.i=l 



N 



ti—to) 



.i=l 



N 

eTT 11/^11- <liminfE^^, 



i=l 



N 



< hm sup Exp 

j->oo 



i=l 
N 



.1=1 



<E. 



AT 



to) 



.1=1 



N 



+ e\{\\h 



i=l 



Since converges to 4> boundedly as jo — > oo, we have 



E, 



TV 



ti—to I 



.1=1 



N 



enil/*ll°o<liminfE 



i=l 



N 



k 

< hm sup "Eixo 



i=l 

N 



.i=l 



<E. 



N 



Umzt,^ 



U=l 



On the other hand, 

- N 

Yi fi{Zti-to] 



E. 



i=l 



E,, 



■ N 

.1=1 



to J 



N 



N 



1=1 



i=l 



From these estimates we conclude that 

- N 

Emu' 



lim E^ 



However, we know that 



lim eJi 



.i=l 



N 



i=l 



^E„ 



:E 



TV 



.1=1 



TV 



.i=l 
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This proves that {Y, P) has the same finite dimensional distributions as 
those for the reflected Brownian motion Z with initial distribution fi. This 
completes the proof of the theorem. □ 

Remark 5.7. By [7], Proposition 3.10.4, Theorems 5.5 and 5.6 hold also 
for step-process approximation defined as 

but with the Skorokhod space D([0, 1],M") in place of the continuous func- 
tion space C([0,1],]R"). 

Now we turn to myopic processes , defined at the beginning of Section 

5. 

Lemma 5.8. Suppose that either (i) /x^ = for every k > 1; or (ii) 
there is a compact subset K of D such that ^j. is a measure on D with 
supfc>i lJ-k{D) < oo and fik{D\K) = for allk> 1. Then the laws o/{X^,P^j 
k>l} are tight in the space C([0, 1],M"). 

Proof, (i) Since, by Lemma 5.2, the laws of {y^, P^^ , ^ > 1} are tight 
in the space C([0, 1],M"') and ^^-fe = X^^_j., we have, in particular, for every 
e > and 5 > 0, that there are /> > and A'^ > 1 so that 

(5.16) P^^. ( sup - > e/i) < 5/2 

for every k> N. 

Recall the oscillation operator osCp[s, t] from the proof of Lemma 5.2 for the 
process X^. As 

oscp[0, 1] < 2 sup osCp[(i - l)2-^ j2-'=] 

iG{l,...,2fc} 

+ sup |Xf2-fe -Xj^s-fel' 

0<ij<2fe : \i-j\2->=<p 

we have 

P^^(osCp[0,l]>e)<P^J sup osCp[(i-l)2-^J2-'^]>e/3) 

MG{i,...,2n / 

(5.17) 

+ PiJ sup |Xf2-.-X^_.|>e/3). 

^0<i,j<2'':\i-j\2-''<p / 
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Since is the reversible measure for the Markov chain {X^^_^.,j = 0, 1, . . . , 2^}, 
for k with < p, 

Pi( sup osCp[(j-l)2-^J2-^]>e/3 
Vie{i,...,2n 

2^ 

<^P^Josc,[(J-l)2-^J2"'=]>e/3) 

= 2^Pj^JosCp[0,2-'=]>e/3) 
<2Sn(L»)Po( sup \Xs-Xt\>e/3 



(5.18) 



which tends to zero as — > oo. Here Pq is the law of Brownian motion X in 
M" starting from the origin. Thus, (5.16)-(5.18) imply that there is A^^i > 
such that 

P^^ (oscp[0, l]>e) <6 for every k>Ni. 

This proves that the laws of {X'^ , P^^ , ^ > 1} are tight in the space C( [0, 1] , M*^ 
(ii) can be established from the tightness of {X^,F^^ , A; > 1} in the space 
C([0, IJjM") by almost the same argument as that for the proof of Lemma 
5.2(ii). So we omit the details here. □ 

Theorem 5.1 is a special case of the following with /x^ = /i = S^xo}^ 
Dirac measure concentrated at xq £ D. 

Theorem 5.9. Suppose that either (i) //^ = ruk for every k>l and /x = 
mo; or (ii) {pk^k >1} is a sequence of measures on D with sup^> /ifc(Z)) < 
oo, Pk{D \ K) = for some compact subset K of D and all k>l, and fi^ 
converge weakly on D to a measure Then the laws of {X^ ,P^^,k > 1} 
converge weakly in the space C([0,1],M") to the reflected Brownian motion 
on D with the initial distribution /x. 

Proof. By Lemma 5.8, the laws of {X^,P^^,A; > 1} are tight in the 
space C([0, 1],M"). Let (Z, P) be any of its subsequential weak limits. Clearly, 
Zq has distribution /i. Let Y be reflected Brownian motion on D with ini- 
tial distribution /x. Denote by Q2 all the dyadic rational numbers. Since 
X^ = for t G [0,1] of the form j2~^, we have from Theorems 5.5 and 
5.6 that {Zt,t G Q2 n [0, 1],P} and {Yt,tG Q2 n [0, 1]} have the same finite 
dimensional distributions. Since both Z and Y are continuous processes, 
(Z, P) must have the same distribution as Y. That is, {Z,P) is the reflected 
Brownian motion on D with the initial distribution //. This proves the the- 
orem. □ 
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